function [output] = fn_C2N(R1,D0i,lambda,y2,s,phi)

rho2 = fn_rho2(s);
rho1 = fn_rho1(s,phi);

output = rho2.*(rho1-lambda*R1)/(1-lambda)*D0i + y2*D0i;

end